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ABSTRACT 

A self-consistent renormalization group flow equation for the scalar A^ 4 theory is 
analyzed and compared with the local potential approximation. The two prescriptions 
coincide in the sharp cutoff limit but differ with a smooth cutoff. The dependence of the 
critical exponent v on the smoothness parameter and the field of expansion is explored. 
An optimization scheme based on the minimum sensitivity principle is employed to ensure 
the most rapid convergence of v with the level of polynomial truncation. 
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I. INTRODUCTION 



In recent years there has been a resurging interest in using the exact renormalization 
group (ERG) formalism for investigating non-perturbative phenomena such as QCD in 
extreme conditions, dynamical chiral symmetry breaking and critical behavior of statistical 
systems. The ERG framework pioneered by Wilson [1], Wegner and Houghton [2], and 
Polchinski [3] was based on the concept of momentum-space blocking transformation 
under which the large momentum degrees of freedom are coarse-grained. Decreasing the 
cutoff systematically then leads to a nonlinear RG evolution equation which dictates how 
the effect of the irrelevant, short-distance modes can be incorporated into the low-energy 
effective blocked action. 

The ERG flow equation may be derived starting from the generating functional for 
the connected Green's functions: 
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where C k = P ^(i-^ P \ P )) * s an additive infrared (IR) cutoff, and Pk,cr(p) is a smearing 
function which approaches <d(p— k) as the smoothness parameter a — > oo. Upon varying the 
effective IR scale k infinitesimally followed by a Legendre transformation, — Wfc[J] + </•$ = 
S k [<&] + ^Q-C^ 1 ®, $ = 8W k /8J, the RG equation for the blocked action S k [3>] is obtained 
as: 
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The above equation has been derived and extensively analyzed [4] [5] [6]. Albeit be- 
ing exact, it cannot be solved analytically and further approximation must be made for 
practical applications. Thus far the most reliable prescription which retains the original 
non-perturbative features has been the derivative expansion [7] where the blocked action 
is written as: 



.(d^y + u k ($) + o(d*) 



(1.3) 



with Zfc($) and U k (<&) being the wave- function renormalization constant and the blocked 
potential, respectively [8]. In the 0(d°) order, the so-called Local Potential Approximation 
(LPA) [9] where Z k (3>) is taken to be unity, the RG evolution is characterized by 
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The equation can now be readily solved by numerical methods [10]. 
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A further approximation to the RG flow is to expand in powers of $ and then 

truncate the series at some order M [11] [12]: 
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thereby transforming the partial differential equation into a set of M coupled ordinary 
equations. Since the field expansion point is arbitrary, it may be chosen at the origin, 
the /c-dependent minimum or any other non- vanishing field values [13]. Thus, polynomial 
truncation not only allows for a direct study of the flow of the coupling constants, numerical 
implementation is also greatly facilitated. However, the drawback of such truncated scheme 
is that it leads to spurious effects such as an oscillation in the value of the critical exponents 
and the appearance of complex eigenvalues at large M [4] [12]. Therefore, it would be 
desirable to employ a prescription which ensures an unambiguous identification of the 
relevant fixed points as well as the convergence of physical quantities. 

In Ref. [14] an optimization procedure based on the principle of minimum sensitivity 
was proposed for measuring the critical exponent v. The parameter a was dialed until it 
reached a value at which the minimum sensitivity condition -^\ a = is satisfied. At the 
optimal a, v depends least sensitively on M due to the maximum suppression of the higher- 
order contributions, leading to a maximum rate of convergence for v. An optimized scheme 
which retains only the relevant set of operators and the leading few irrelevant operators to 
characterize the RG flow can have non-trivial consequences in more complicated systems 
such as fermionic or gauge theories. In particular, the computation could be simplified 
enormously if certain classes of higher-order Feynman diagrams are suppressed. 

In the present work we continue our search for an optimized RG prescription at the 
level of LPA. However, instead of following the conventional ERG approach as previously 
considered [14], we start from the one- loop contribution and introduce a "multiplicative" 
regulating smearing function Pk,a(p) in the momentum space integration: 



Ui 1] m = ~ 2 jhAp) ^ [l + • (1-6) 



v 



The use of pk,a(p) allows one to see clearly how each Feynman diagram is regularized. From 
the above equation, a self-consistent RG equation can be constructed by a differentiation 
with respect to k followed by the substitution V"(&) — > U%($>) on the right- hand-side: 



p2 



(1.7) 



This procedure is analogous to the Schwinger-Dyson self-consistent improvement. The 
main focus of this paper is to illustrate how Eq. (1.7) can provide an alternative, viable 
prescription that retains the non-perturbative information contained in the original full 
RG equation for iSfc[<&]. In particular, the convergence of physical quantities such as the 
critical exponents extracted using polynomial truncation will be demonstrated. 
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The organization of the paper is as follows. In Sec. II we briefly review the essen- 
tial features of LPA and compare with the self-consistent RG equation, Eq. (1.7). The 
sharp cutoff limit of the RG flow is discussed in Sec. Ill and the dependence of the crit- 
ical exponent v on the field expansion point $o an d the level of truncation M is also 
explored. We expand the potential about its fc-dependent minimum, a prescription which 
does not respect the Z 2 symmetry of the original flow equation. Nevertheless, the resulting 
v converges remarkably with increasing M up to M = 26. In Sec. IV the RG equations 
generated by smooth smearing functions are presented. Unlike the sharp cutoff case, sin- 
gularities which resemble spontaneous symmetry breaking are present in the self-consistent 
integro-differential equation for the fixed-point solution U*($>) when employing a smooth 
cutoff. To overcome the complication we again make use of the non-invariant expansion 
and obtain numerical solutions for v up to M = 20. Sec. V is reserved for summary 
and discussions. In Appendix A we summarize some detail associated with polynomial 
truncation of the sharp flow equation. The manner in which propagators are modified by 
LPA and the self-consistent RG schemes is discussed in Appendix B. In Appendix C we 
explore the large a limit of the self-consistent RG. 



II. LPA AND SELF-CONSISTENT RG 



We first give a brief recapitulation of the basic features of LPA. The ERG equation 
can be constructed by introducing a smearing function pk,a(x) which governs the coarse- 
graining procedure and k acts as an effective IR cutoff, separating the low- and the high- 
momentum modes. The propagator is modified as: 

A( P ) = i _ AkAp) = i^M = fccfE), (2.1) 

where Pk,a{p) + Pk,a(p) = 1- In LPA where Zfc(<&) = 1 the evolution is characterized by 
the flow of C/fc(<]>). For a generic smooth function Pk,a(p) = ®a(k,p), it is convenient to 
write Q a (k,p) = t where < t < 1, and the reparameterization leads to 

dUm _ Sa f^=^ p d (t) U»{*) 

dk ~ 2 y t(p=0) P i(t) + (i-t)u^y {2 - 2) 

Thus, by inverting Pk,a{p) for p(t), the RG equation of is readily obtained. 

Clearly, the evolution equation depends on the choice of the smearing function. While 
a sharp cutoff provides a well-defined boundary between the high (integrated) and the 
low (unintegrated) momentum modes and yields a partial differential equation, no clear 
separation exists for a smooth cutoff and the RG equation retains its integral-differential 
form. In principle, physical content extracted from solving the RG flow should not depend 
on the shape of the cutoff which merely reflects how the irrelevant degrees of freedom 
are integrated over. However, due to the approximation employed such dependence is 
explicitly present and the critical exponent v exhibits a small dependence on a [7]. The 
dependence may indicate that the effect of certain relevant operators has been neglected 
[15]. 
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Within the framework of LPA, we see that Eq. (1.4) may also be obtained by a similar 
self-consistent "trick" starting from [16]: 



In 



(2.3) 



This can be understood by noting that since we are now using an "exact" propagator, all 
the higher-order Feynman graphs that belong to the one- loop resummed class (e.g., daisy 
and super-daisy) are automatically incorporated. More generally, we may write: 
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(2.4) 



(2.5) 



Here Vk,a{p) acts as a projection operator on the massless propagator 1/p 2 and turns it 
into an effective propagator pk,a{p)/p 2 which contains only the high momentum modes 
p > k. For example, in the sharp limit one would have A^(p) = @(p — k)/p 2 . 

Let us examine how the alternative, self-consistent RG prescription based on Eq. 
(1.6) differ from the usual LPA. First we note that since the coarse-graining procedure is 
initiated in the high-momentum regime where V"(Q)/p 2 <C 1, Eq. (2.3) may be expanded 
as 
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Differentiating both sides with respect to In k and invoking the self-consistent prescription 
then leads to 
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(2.7) 



In the first term of the equation above one finds a complete factorization between the 
running of the scale k and the $ dependence. This is the equation which we shall focus 
on in the subsequent sections. The second term, representing the difference between the 
self-consistent RG and the LPA, can also be written as: 
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From Eq. (2.6) or (2.8), we see that the two RG schemes coincide when Pk,a(p) ^ 1 or 0. 
That is, AU k ($) = when V% jtT (p) = V k , a (p) or p kj<J (p) = pk,*(p)- However, this holds 
only for the trivial case of a sharp cutoff; smooth cutoffs such as p~k,b(p) = 1 — e~ al <Pl k ) 
and Pk, m (p) = (p/k) m / (1 + (p/k) m ), on the other hand, give a small deviation for large m 
and b in the regime p/k ^> 1. 

The salient feature of our approach is that once the perturbative results are known, 
they can be readily transformed into the corresponding self-consistent RG equations. Gen- 
eralization to higher order derivative expansion can be achieved with ease. For example, 
at order 0(<9 2 ), we have: 
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and [17] [18] 
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(2.10) 



We comment that the above prescription, being reminiscent to the Schwinger-Dyson ap- 
proach, has been widely employed to treat more complicated issues such as the effective 
quark propagator [19] and dynamical chiral symmetry breaking in gauge theories [20]. 
One great advantage here is that, unlike the exact RG approach, the flow for Z k (&) is 
completely analytic, even in the sharp cutoff limit, and the numerical algorithm can be 
readily implemented. Incorporating the effect of Z k (<&) with the above approach has been 
shown to improve the accuracy of v [21]. In what follows we shall set Z k (3>) = 1 (77 = 0) 
for simplicity. 

Although we shall mainly be concerned with smearing functions that approach a step 
function as a — > 00 in our self-consistent RG scheme, other type of regularization can also 
be utilized. For example, one could choose the smearing function to be of the Pauli-Villars 
type: 

~ , ^ / A 2 \ d/2 ( k 2 



In a similar manner, the spirit of the dimensional regularization can also be encapsulated 
by using p e (p) = d~ e p/(2n)~ e in integer d; with A — > Xp, e , the RG flow is described by p, 
the renormalization point. 



III. SELF-CONSISTENT RG - SHARP CUTOFF 

Before considering smooth cutoffs, it is instructive to first review the sharp cutoff 
results [14]. In order to recover this limit, caution must be taken to invert t = 0(1 — z). 
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Here it must be interpreted that t changes continuously from 1 to at z = 1, thereby 
leading to [4] : 



dt In 



1 + 



= -In 



l+uzm , (3.i) 



where U k ($) = C 2 &~ d t4($), $ = Cfc~ (d ~ 2)/2 $, with C = VWKi = ^ (An) d / 2 T(d/2). The 
fixed-point solution is given by 



(d - 2)$U*' + dU* = -In [1 + U* 



(3.2) 



with U* (0) = from reflection symmetry. 

However, it has been pointed out that only certain values of = U*(0) = — ln(l + 
gfe)ld can lead to physically meaningful solutions [4] [11]; in all other cases a singularity 
of the form U* ~ ln($ c — $) will be encountered and the potential becomes complex 
beyond some critical value $ c . In d = 3, only two solutions have been found; the first 
one is = which corresponds to the trivial Gaussian fixed point and the second one, 
% = 0.206343 = s (or g\ = -0.461533 = s), the Wilson-Fisher. The qualitative feature 
of the fixed-point solutions is illustrated in Fig. 1. 
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Figure 1. Singularity <J> C as a function of go in d=3. Notice that there are only two physical solutions, 
one for the Gaussian and the other one for the Wilson-Fisher fixed point. 



As we mentioned in the Introduction, instead of solving the full flow equation, we 
may expand £/&($) in powers of $ and terminate the series at i> 2M , thereby turning 
the nonlinear partial differential equation into a set of M coupled ordinary differential 
equations. We consider below expansion schemes that are both Z 2 -invariant and non- 
invariant. 



a. invariant expansion 
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The simplest invariant polynomial expansion of E/fc(<&) can be carried out around 
$ = 0, as in Eq. (1.5). The (3 functions, n = kdg n /dk, can be written as: 



02 

h 

06 



-2g 2 - G 4 , 

-m - (-3Gl + G 6 ). 



(3.3) 



(2d - 6)g 6 - (30G| - 15G 4 G 6 + G 8 ) : 



where G n = g n /(l + #2) [22], and the fixed points can now be located by setting n = 
for all n. If only g 2 and g 4 are kept, we find, besides the trivial fixed point (g^dt) = 
(0,0), an additional Wilson-Fisher fixed point (g^gl) = (—q^, (6+e) a ) • Linearizing the 
flow about the above fixed point, the eigenvalues can easily be found to be [(—3 + 2e) ± 
V9 + 6e + 7e 2 ]/3, or (—1 ± v^22)/3 for e = 1. By our convention, the largest negative 
eigenvalue corresponds to the inverse of the critical exponent Thus, we have v « 

0.5272 at M = 2. We remark that in the spirit of the conventional e-expansion where 
e is first taken to be "small," the matrix element d0 4 /dg2 which is 0(e 2 ) is discarded, 
and v = 0.6 instead [23]. Thus, when characterizing the flow by powers of there is a 
reshuffling of contributions from the e series. 
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Figure 2. Oscillation of critical exponent v and g% as a function of M. 



However, it is well-known that v obtained via polynomial truncation with a sharp 
cutoff generally does not converge but oscillates with M, as depicted in Fig. 2. We find 
that at M — 22 an additional relevant direction is generated; and as M is further increased, 
more spurious solutions are found and the identification of the Wilson-Fisher fixed point 
becomes ambiguous. The oscillation in v is intimately linked to the oscillation in g\ and 
the cause is related to the singular structure of the non-truncated fixed point potential 
[/*($) on the complex plane [4]; namely, the asymptotic behavior is dominated by a pole 
$ c = |$ c |e l6lc having 9 C « ±7r/2 that gives rise to an approximate four- fold periodicity. 
Such truncation dependence is not unexpected since the fixed point(s) one identifies at 
each order is only approximate and not the true Wilson-Fisher fixed point; therefore, the 
would-be irrelevant set of operators will continue to evolve, leading to the M dependence 
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of physical quantities. In order to improve the convergence of v in polynomial truncation, 
Aoki et al has carried out an expansion about a ^-invariant field xo = ^o/2 [13]: 



M 



n=0 



^(X) = E (X - XoT, c n (k) = Ujr\ x0 ) = ^ _ (3.4) 



X=XO 



The convergence has been found to improve dramatically. A more thorough discussion of 
this expansion scheme is presented in Appendix A. 

b. non-invariant expansion 

Alternatively, one may expand the potential around some non- vanishing 

&(*) = E*5T(*-*T. W'l^'^l^u/ < 3 ' 6 > 



n=0 



Here the presence of odd powers of $ will explicitly break the Z 2 symmetry. However, this 
novel prescription becomes particularly useful when spontaneous symmetry breaking takes 
place, or when the running mass parameter g 2 (k) is still negative during the RG evolution. 
It also allows us to monitor the flow around $' instead of the origin. 

Choosing = i> , the /c-dependent minimum and using Eq. (3.5), we obtain 



(3.6) 



$o = -(l-|)$o + G 3 /g 2 , 
h = 93^>o - 2g 2 + (1 - |)$o#3 + (Gl - G 4 ), 

= 9^0 - (1 + ^93 + (1 - l)$o9± + (-2GI + 3G 3 G 4 - G 5 ), 
fa = g^ Q - eg 4 + (1 - |)$ ^5 + (6<?£ - 12G 2 3 G 4 + ZG\ + 4G 3 G 5 - G 6 ), ■ ■ ■ 

where G n = At the fixed point where i>o vanishes, we have $q = ( 2 -e)^*(i+g*) • 

When solving the coupled equations j3 n = to locate the fixed points, one finds that the 
number of spurious solutions again increases with M. Moreover, due to the loss of Z 2 
symmetry, the Wilson-Fisher fixed point no longer has just one relevant eigen-direction, 
as would be in the ^-symmetric case. 

The new Wilson-Fisher fixed point which dictates the critical behavior of the Ising 
universality class satisfies the following conditions: (i) all the fixed-point coupling constants 
<7*'s are real, (ii) the extremum of i7*(i>) is close to $q = 1-9287, the field value obtained 
from solving the full non-truncated LPA equation, and (iii) at M = 2j (j > 2) there are 
£ = j + 1 relevant eigen-directions which, by our convention, have negative real part for the 
eigenvalues. Let us see how the above criteria work. First we note that g 2 < for M < 6, 
an indication that $o is a local maximum. At M = 4 there are only three non-trivial 
solutions (g 2 ,g 3 ,g 4 ) = (-£,0, g_) and (-0.4125, ±0.2327, 0.4456), with the first being the 
usual ^-symmetric solution at $o = 0. We choose the asymmetric solution with g 3 < 0, 
which yields $0 = 1-9203 and v = 0.4229. For M > 6, the extremum becomes a local 
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minimum and we have I = j + 1. This is due to the fact that the fields that break the Z2 
symmetry remain relevant throughout. In addition, there are two more coupling constants, 
g\ and <?|, which are relevant by power counting. For example, at M = 10 the relevant 
set should be {92,93,94,951 97, 99}, i-e., I = 6. Note that the classification of operators 
around the trivial Gaussian fixed point is still based on the counting of their canonical 
dimensions [g n ] = d — n(| — 1). 

The fixed-point potential £/*($) is plotted in Fig. 3 as a function of $ for various M. 
Due to the Z2 asymmetry and the choice $0 > as the expansion point, only £7*($) with 
positive $ is tracked. We note that the approximation first improves with increasing M 
and nearly coincides with the full non-truncated potential at M = 20. Incidentally this 
is the order at which the term linear in $ that represents the leading symmetry-breaking 
contribution is minimized. However, beyond M — 20 the potential begins to deviate away 
from the non-truncated LPA solution. At M = 28 a sensible result can be obtained only by 
choosing a complex eigenvalue. The remarkable convergence of v with increasing M up to 
M = 26 is depicted in Fig. 4. Note that in the Z2-invariant polynomial truncation scheme, 
depending on the point of expansion, v is found to converge to or oscillate about 0.689, 
which is the value obtained by solving the non-truncated LPA equation. On the other 
hand, our non-invariant prescription gives 0.649 for a large range of M. Such discrepancy 
is not unexpected since the symmetry of fixed points for the Ising universality class differs 
in the two cases. It is interesting also to note that the result is very close to the value 0.651 
obtained via the optimization procedure in which the most optimal smoothness parameter 
is sought to achieve minimum sensitivity [14]. 
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Figure 3. Fixed-point potential U* ($) obtained by polynomial truncation. Exact LPA solution is 
shown as solid black line. The agreement with the full LPA solution first improves with 
increasing M and nearly coincides at M =20. 
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Figure 4. Critical exponent v as a function of M obtained by an non-invariant expansion of U* (<£) 
about the minimum &q. 

We comment that besides the origin or the /c-dependent minimum $o, i n principle $' 
can be arbitrarily chosen. However, the rate of convergence is closed linked to the choice 
of Aoki et al has shown that expanding the sharp cutoff RG equation for £?*($) about 
the minimum, $0, yields a larger radius of convergence compared to at the origin [13]. 
In fact, it leads to the exact solution in the large N limit. Thus, when expanding about 
the minimum $o, a sharp cutoff smearing function with a —>■ oo would be the optimal 
"smoothness" that satisfies the minimum sensitivity principle, as we shall see in the next 
Section. 



IV. SELF-CONSISTENT RG - SMOOTH CUTOFFS 

In this Section we examine the fixed-point potential £/*($) obtained via polynomial 
truncation with smooth cutoffs. Unlike the sharp cutoff case, the RG flow is characterized 
by the following (dimensionless) integro-differential equation: 



k^---(d-2)$^ + d 
dk 2 V ' ^ 



d_ 



In 



z 2 + uzm 



(4.1) 



- / dt z d (t) In z 2 (t) 
Jo 



The smearing functions are chosen in such a way that the conditions t(z = 0) = 1 and 
t(z = oo) = are always satisfied. However, we shall find that within a certain range 
of k in the smooth case the running mass parameter §2(k) (as well as g%) are negative 
quantities and singularity is encountered in j3 n as the denominator z 2 {t) + &£($) goes to 
zero. 
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The reason why such singularity is absent in the sharp cutoff but present in the 
smooth cutoff for self-consistent RG can be understood from the manner in which the 
large momentum modes are eliminated. In the case of exact RG a sharp cutoff provides 
a clear separation between the integrated (z = p/k > 1) and the unintegrated modes 
(z = p/k < 1). With a negative mass curvature g\ = s ~ —0.461533 (in unit of k) for 
the non-truncated potential, the integration is completely analytic because the condition 
1 + s > is always satisfied in the denominator. However, in the self-consistent RG 
formulation no sharp boundary exists for a smooth cutoff and for any finite a the range 
of the z integration is stretched from zero to infinity. Therefore, singularity sets in when 
z 2 + s < 0. The situation confronted here actually resembles that of spontaneous symmetry 
breaking where the potential is also double- welled. However, we emphasize that £/"*($) 
describes a scale-invariant theory in the symmetric phase; as we integrate down to k = 
in the dimensionful RG equation, the potential will become single-welled with a unique 
minimum at the origin [10]. True symmetry breaking is marked by lim/^o 92(k) < 0. 
Notice that the integration with a smooth cutoff would be completely analytic had the LPA 
prescription been employed; the difference between the two schemes may be attributed to 
the manner in which the propagators are modified. The details are further discussed in 
Appendix B. 

To overcome the difficulty associated with the divergence near z = \f—^fa, we expand 
the potential about a non-zero minimum $o where U"($>o) = 92 > 0, as in Eq. (3.5). In 
analogy to Eq. (3.6), we have, for the smooth case: 

io = -(l-|)$o + f z d (t)g 3 /g 2 , 

e - f 1 

fo = gs&o - 2g 2 + (1 - -)$o£ 3 + J dt z d {t) (0f - 

h = gSv - (i + \)gz + (i - |)$o£ 4 + ^ dt z d {t) (-2<?f + 3^4 - 

h = h®o - m + (1 - ^)$o#5 + f dt z d (t) (6Gl - YlQlQ^ + 3Gj + 4Q 3 Q 5 - <? 6 ) , • • • 

z Jo 

(4-2) 

where Q n = zi ffi + § 2 ■ The absence of linear term implies $o = (2-e)g 2 Jo dtz d (t)Q 3 . The 
flow of $o is connected with that of g n with n odd; naturally all odd coupling constants 
vanish when $o = 0. 

Consider the two smooth representations: 

p k M=e-^ > ^4t)=(—) 1/ \ (4.3) 

and 

/I \ l/m 

Pk, m (p) = [i + (p/krr 1 — z{t) = (- - 1) . (4.4) 

We choose a = In 2 so that Pk,b{p = k) = 1/2 and the smoothness of the cutoff is con- 
trolled by a single parameter b [24] . From Figures 5 and 6 depicted below, it is clear that 



12 



the sharper the cutoff the more rapid the convergence for v in the polynomial truncation 
scheme. That is, m, b = oo is the optimal value at which the minimum sensitivity condition 
jM\b,m — is met. As b(m) is decreased, v oscillates with increasing amplitude. Com- 
paring the two figures, one also notices a smaller oscillation amplitude for the exponential 
cutoff. 



0.8 



0.75 
0.7 



0.65 



0.6 

6 8 10 12 14 16 18 20 

M 

Figure 5. Critical exponent v as a function of b for various M. The optimal smearing function 
corresponds to a sharp cutoff with b=oo. As b is decreased, v begins to oscillate with 
increasing amplitude. 
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0.75 
0.7 



0.65 



0.6 

6 8 10 12 14 16 18 20 

M 

Figure 6. Critical exponent v as a function of m for various M. The "optimal" smoothness here 
again is m=oo. 



We emphasize again that the optimal smoothness a obtained from the minimum sen- 
sitivity principle generally depends on the model under study as well as the approximation 
employed. While a sharp cutoff is the optimal smoothness when $' is chosen to be the min- 
imum $o °f the potential, for an arbitrary the optimal a which satisfies the minimum 






m = 1000 




m = 30 
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sensitivity condition would correspond to a smooth cutoff. On the other hand, a sharp cut- 
off yields a rather poor convergence for $' = 0. A general criterion of finding the optimal 
a can be found by analyzing the massless blocked propagator A C7 (z) = 1/P%(z) = f>k,a/z 2 , 
as has been recently shown by Litim (see Appendix B) [15]. 

V. SUMMARY AND DISCUSSIONS 

In the present work we introduced a self-consistent RG prescription starting from the 
one-loop perturbative expression. Our self-consistent equation coincides with LPA in the 
sharp cutoff limit but differs with a smooth cutoff. The advantage of this novel approach 
is that the perturbative independent-mode results can be readily "dressed" into a set of 
coupled non-linear RG equations. Complication arising from taking the sharp cutoff limit 
in the derivative expansion can also be avoided completely. 

Due to the absence of a sharp boundary in the momentum integration in the smooth 
cases, a singularity is encountered in solving the truncated self-consistent equations, as the 
inverse propagator z 2 + g\ vanishes. We circumvent this difficulty by expanding U*(&) 
about its minimum $ where g\ > 0. Such expansion is non-invariant since it contains odd 
powers of $ that destroys Z^ symmetry. The system remains, however, in the symmetric 
phase; as we integrate down to k = in the dimensionful flow equation, the moving 
minimum $>o(k) vanishes. Nevertheless, the formal equivalence provides a useful tool for 
exploring systems in the state of broken symmetry. 

In the absence of Z^ symmetry, the physical, Wilson-Fisher fixed point which governs 
the critical behavior of the Ising universality class may be readily identified from the 
following conditions: (i) all the coupling constants g*'s are real; (ii) the potential has 
an extremum close to $q = 1.9287, the value obtained without truncation; and (iii) the 
number of relevant eigen-directions is j + 1 for M = 2j (j > 2). Having identified the 
fixed point, the exponent v is extracted from the truncated polynomial^ series. As shown 
in Sec. Ill, when a sharp cutoff is employed, v oscillates with M when [/*($) is expanded 
about $ = 0. On the other hand, the convergence improves dramatically in the non- 
invariant expansion about the minimum <E>o- We truncated the series at M = 28 where 
the eigenvalues become complex. From Figure 4, we see that with increasing M the series 
converges rapidly to approximately v = 0.649, which differs from the full LPA solution 
0.689. The discrepancy may be attributed to the difference in the symmetry of the two 
approaches. In Sec. IV where smooth cutoffs are considered, one finds that the amplitude 
of oscillation grows with a smaller a, and the rate of convergence decreases. Thus, the 
sharp cutoff is the "optimal" smearing function for the self-consistent RG expanded about 
the minimum $q- Indeed, the minimal sensitivity condition -^\ a = is reached as a — > oo. 

The results obtained in this paper demonstrates that when solving the truncated RG 
equations for the convergence of physical quantities is strongly influenced by two 

factors: (i) the expansion point <&q, an d (ii) the smoothness parameter a in the smearing 
function Pk,a(p)- An optimal choice of $o? usually the minimum of Uk(Q), yields the 
largest radius of convergence in the expansion [13]. However, when $ is n °t optimal, one 
can fine-tune a so that physical quantities such as v would exhibit minimal sensitivity on 
the order of truncation M. 

The optimization prescription presented here can be readily extended to more com- 
plicated systems. Polynomial truncation applied to models with fermionic interaction can 
be considered [25]. A self-consistent RG prescription with a smooth cutoff can also be 
formulated via the operator cutoff regularization in the proper-time space [21] [26]. In 
order to test the success of our self-consistent RG scheme, the effect of wave function 
renormalization constant Zk(&) must be incorporated. The results will be reported in the 
forthcoming publication [27]. 
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APPENDIX A: SHARP RG FLOW BY POLYNOMIAL TRUNCATION 

In this Appendix, we provide some detail of the RG behavior near the Wilson-Fisher 
fixed point, based on the polynomial truncation. The evolution equation obtained with a 
sharp cutoff reads, in d = 3: 



There are several ways to obtain an approximate solution to the fixed-point equation. 
We summarize two schemes below: 

1. expansion about i> = 0: 

Here we assume the potential to be of the form: 



Brazil. 



- -®u' k m + = -in [i + u£(*j . 



(A.l) 




M _ 



d 2n U k 




(A.2) 



and the flow is characterized by a set of M coupled equations for g 2 



(1) M = 2: 



At this order, we have 



(3 2 = -2g 2 - 



1 + 02 ' 
3gi 



(A.3) 



At = ~94 + 



(1 + 92) 2 



+ 0(g 6 ), 



and the Wilson-Fisher fixed point is located at (g^dt) = ( — \, !§)• Linearizing the flow 
then yields 




(A.4) 
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where 



M = 



002 


d§2 










dg 2 


dg4 



-2 + 



9a 



(l+3 2 *) 2 

(T+IJF 



-e + 



1+32 



-2 + f 



4ei 
6+e 



6+e 
6 



(A.5) 



is the linearized RG matrix about the Wilson-Fisher fixed point. Upon solving the eigen- 
value problem for M., we find v = 0.527 which deviates from the exact value 0.689, found 
by solving the non-truncated equation (A.l). The discrepancy is due to the fact that 
our approximate Wilson-Fisher fixed point still lies too far away from the exact solu- 
tion (s 2 , S4) = (—0.4615, 0.4970). Thus in the neighborhood of our approximate solution, 
higher-order operators continue to evolve and therefore must be included in order to im- 
prove the accuracy. 

(2) M = 3: 

Taking into consideration the running of c/q, we have 

#4 



P2 

h 
fa 



-2g 2 - 
-94 + 



I + 92' 



96 



(l + g 2 ) 2 I + 22' 

+ 0(g 8 ), 



(A.6) 



30^2 15^6 



+ 



(l + <72) 3 (I + S2) 2 



which leads to 



(— * — * — *\ 
92, 94, 96) 



1 4 16 
v 3' 9' 27^' 



(A.7) 



A comparison with the M = 2 calculation shows that the fixed points are now closer to the 
exact solution with the inclusion of g~Q, and hence a more accurate v is expected. Indeed 
for M = 3 we find v = 0.585. The complex singularities near the origin at this order are 
found to be $ c = ±0.58998 ± 2.20183i, or (|$ c |, 6 C ) = (2.2795, ±0.4167tt). 

As discussed in Sec. Ill, the terms from the e expansion are reshuffled in the polyno- 
mial truncation scheme. Retaining terms up to 0(e 3 ) at this order, we find 



M = 



f-2 



\ 



+ i + 

2e 2 



5£ 
9 



+ 

14e 3 



55e 3 
54 



-(i + 

e + 



e. _i_ be 
6 18 
10e 2 



55£ 
108 



10e d 



3 

10e 2 



+ 



55e d 



) N 

K 1 ~ 6 ~ 18 T 108 / 



70e d 



2 + 3e + 



25^ 



+ 



108 
275e 3 
18 



(A.8) 



which, for e = 1, leads to 



v = 



or v m 0.653. Comparing Eq. (A. 9) with that obtained in 



(A.9) 



(A.10) 
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we see that they agree up to 0(e); the discrepancy in the e 2 term is due to the negligence 
of certain Feynman diagrams as well as the wavefunction renormalization. 

(3) M > 4: 

With increasing M, our approximate fixed point moves even closer to the non- 
truncated solution and a greater accuracy is attained. The results are summarized as 
follows: 



M 


#2 


— * 

94 


— * 

9 6 


#8* 


— * 

#10 


#1*2 


#1*4 


— * 

#16 


2 


-.1429 


.2449 














3 


-.3333 


.4444 


.5925 












4 


-.4664 


.4977 


1.1273 


2.7804 










5 


-.5163 


.4995 


1.3056 


4.2462 


9.9214 








6 


-.4508 


.4952 


1.0675 


2.3609 


-2.303 


-77.3326 






7 


-.4349 


.4915 


1.0048 


1.9537 


-4.2865 


-84.4919 


-187.96 




8 


-.4665 


.49776 


1.1278 


2.7837 


0.0191 


-66.3291 


248.552 


18590.5 



Table 1. j* as a function of M. 

The relation between different #*'s can be obtained by expanding U* about $ = 0: 

(A.ll) 



oo 



#2r 



0-(») = -!ln(l + jS ) + f *»+ , 



2n 



For large M, the value of g\ approaches s = —0.461533 and g^n can be rewritten in terms 
of s as: 

gt = -2s(l + s), 
# 6 *=2s(l + s)(l + 7s), 

g* = _60s 2 (l + s)(l + 3s), (A.12) 
#* = 40s 2 (l + s)(2 + 43s + 83s 2 ), 
g* 12 = -80s 2 (l + s)(-2 + 45s + 549s 2 + 880s 3 ), • • • 

From the non-truncated RG equation we see that the potential becomes singular when 
the argument inside logarithm vanishes, i.e., 1 + U* ($) = 0. Therefore, when expanding 
ln[l + U* ($)] in power series of for large enough M and g\ ~ s, one observes a four- 
fold periodic oscillation of (?* in (+, +, — , — ) pattern, due to the presence of a complex 
singularity near the origin whose phase is approximately tt/2. Thus, increasing M will at 
first improve the result by pushing the singularity further away from the origin $ = where 
the expansion is made, but eventually the series starts to oscillates about v = 0.689 ±0.008. 



2. expansion about $ 2 7^ 0: 



17 



In order to improve the convergence in polynomial truncation, Aoki et al has consid- 
ered an expansion about some non- vanishing field instead of the origin [13]. Employing 
X = & 2 /2 which preserves the Z 2 symmetry of the original action, Eq. (3.1) becomes 

l + ^(x) + 2x^(x)]- (A.13) 
The expansion 



k^-(d-2)x— + d 



dk 



Uk(x) = "In 



M 



u k (x) = J2 c -^(x-x\k)r, 



Cn(k) = Ui n \ X ') = 



d n U, 



n=0 



dx r 



(A.14) 



leads to 



Co = X'ci - (4 - e)c + (2 - e) X 'ci - ln(l + c x + 2 X 'c 2 ), 
ci = x'c 2 - 2ci - x' [(-2 + e)c 2 + 2c 3 ] - 3c 2 , 

c 2 = X'c 3 - (4 - e)c 2 + (2 - e)(2c 2 + X 'c 3 ) + (3c 2 + 2 X 'c 3 ) 2 - (5c 3 + 2 X / c 4 ), 
c 3 = X'c 4 - (4 - e)c 3 + (2 - e)(3c 3 + X 'c 4 ) ~ 2(3c 2 + 2 X / c 3 ) 3 , 

+ 3(3c 2 + 2x'c 3 )(5c 3 + 2 X 'c 4 ) - 7c 4 + 0(c 5 ), 
c 4 = (4 - 3e)c 4 + 6(3c 2 + 2 X 'c 3 ) 4 + 28(3c 2 + 2 X / c 3 )c 4 - 12(3c 2 + 2 X / c 3 ) 2 (5c 3 + 2 X 'c 4 ), 

+ 3(5c 3 + 2 X / c 4 ) 2 +0( C 5 , C 6 ),--- 

(A.15) 

where c n = c n /(l + c\ + 2%'c 2 ). As illustrated in Sec. Ill, a more rapid convergence 
would be obtained if the expansion is made around ore near the minimum xo = 1-844. We 
summarize in Table 2 below for the results obtained about x' = 2: From the Table we see 
a remarkable convergence of v toward the full LPA solution in this expansion scheme. 



M 


V 


Cl 


C2 


C3 


C 4 


C5 


C6 


C7 


C8 


C9 


3 


.66781 


-.0084 


.2664 


.0829 














4 


.68698 


.0303 


.3343 


.1091 


.0137 












5 


.68803 


.0490 


.3586 


.1155 


.0145 


-.0071 










6 


.68907 


.0512 


.3619 


.1167 


.0151 


-.0072 


.0016 








7 


.68956 


.0498 


.3601 


.1162 


.0150 


-.0066 


.0016 


.0026 






8 


.68952 


.0492 


.3591 


.1159 


.0148 


-.0067 


.0011 


.0023 


-.0040 




9 


.68948 


.0491 


.3590 


.1159 


.0148 


-.0066 


.0011 


.0025 


-.0040 


.0012 



Table 2. The critical exponent v coupling constants c„ as a function of M. 



18 



0.7 



0.69 - 



^ 0.68- / 
0.67 / 

0.66 ' 1 1 1 1 1 

3 4 5 6 7 8 9 

M 

Figure 7. Critical exponent v as a function of M obtained using the invariant expansion about x'=2. 

APPENDIX B: MODIFIED BLOCKED PROPAGATORS 

We have seen that with a smooth cutoff the momentum integrations for the fixed-point 
solution are divergent in the self-consistent RG formalism, while they remain completely 
analytic in the LPA. The distinction can be readily understood from the manner in which 
the propagators are modified in the two approaches. 



a. exact RG 



In the exact RG approach, the second functional determinant is given by S'l\Q\ = 
C^ 1 +S"[$>] where C^ 1 = p 2 (l— Pk,a) / Pk,a- This in turn implies a modified (dimensionless) 
propagator of the form 

where z = p/k, P 2 (z) = z 2 / ' pk,a and #2 = U'^(0). Let us first examine the behavior of 
P 2 (z). For definiteness we take the smooth smearing functions to be pb{z) = 1 — e~ az as 
considered in Sec. IV. 

We see that as z is decreased, with the exception of b = 2, P 2 attains a minimum and 
then diverges as z — > 0. The minimum zq, located by dP 2 /dz = 0, satisfies the relation 
ZQe~ az "/(l — e~ az ") = 2/ab and implies a minimum P 2 (zq) = 2z 2 ~ b e az o jab. As can be 
seen from Table 3, P 2 (zq) approaches unity in the limit b — > oo. As recently addressed by 
Litim [15], b can be chosen in such a way as to achieve the maximum radius of convergence 
R = lim n ^oo a n /a n +2 in the amplitude expansion whose coefficients are given by 



Jo 



dzz d+1 (k^^)p- n (z). (B.2) 
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By using the following criterion [15]: 



.Ropt = max^min z > P%{z)j , (B.3) 

the optimal smoothness parameter associated with this particular class of smearing func- 
tions is given by 6 opt = 2/ln2, with R opt = 2. On the other hand, based on the principle of 
minimum sensitivity, it was found that b* = 3 yields the fastest rate of convergence with 
respct to polynomial expansion for v [14] . We emphasize that the two optimization proce- 
dures are inherently different; while 6 op t is a universal result for the exponential smearing 
function, the vale of b* depends on the model and the physical quantity under study [28]. 



b 


zo 


P b 2 (z ) 


2 





l/ln2 = 1.4427 


2/ln2 


1 


2 


3 


1.0324 


1.9974 


4 


1.1603 


1.8821 


5 


1.1849 


1.7508 


20 


1.0861 


1.2122 



Table 3. Location of the minimum of P£ as a function of b. 

One may also consider the power-law cutoff p m {z) = z m /(l + z m ). From Table 4 we 
may readily conclude that m op t = 4. The minimum sensitivity condition utilized in Ref. 
[14] gives m* = 5.5. 



m 


zo = (^)V- 


D2 ( „ \ _ (m-2\2/m( m ^ 
^m(. z 0) ~ { 2 > y m -2) 


2 





1 


3 


.7937 


1.8899 


4 


1 


2 


5.5 


1.1071 


1.9261 


20 


1.1161 


1.3841 



Table 4. Location of the minimum of as a function of m. 

In Fig. 8 the blocked propagator is depicted for the exponential cutoff. Note that for 
6 = 2, 1/P£(z) acquires an effective mass which remains finite in the limit z — > 0. 
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Figure 8. 1/P^(z) as a function of z with a=ln2 and b=2,3,4,10,20. The propagators is IR finite as 
z^O. 



Let us now consider the effect of g 2 . In the symmetric phase where g 2 > 0, the change 
is only a shift P 2 — > P 2 + g 2 . On the other hand, if g 2 < as in the case of spontaneous 
symmetry breaking or during the course of RG evolution where g 2 is still negative, the 
modified propagator can become singular if P b 2 (z) + g 2 < 0. However, this does not happen 
in LPA where —g 2 = 0.461533 < minP 2 (z) = 1. Thus in the LPA prescription with a 
smooth cutoff, the momentum integrations are completely analytic. We comment that in 
the case of spontaneous symmetry breaking where g 2 remains negative down to k = 0, one 
expects the inverse (dimensionful) propagator (p 2 /pb) +#2 to vanish at some scale fc cr , and 
the region below k < k cr becomes inaccessible by conventional perturbation theory 



b. self-consistent RG 



Let us now turn to the self-consistent RG prescription. In this formalism, the modified 
propagator reads 

A a (z) = 44^- (B.4) 

Z +92 

Contrary to Eq. (B.l), §2 is also influenced by the presence of p a (z) and must not be 
treated separately. Using the exponential smearing function pb(z) as an example, we find 
that for <j2 > the qualitative feature of the modified propagator remains the same except 
for b — 2; while P b 2 approaches l/ln2 as z — > in the exact RG formalism, it diverges in 
the self-consistent RG prescription. The minimum of P b (z) is located by solving 



zi + g 2 _ 2zl~ b e az o 



1 _ e ~az° a b 



(B.5) 



The major difference between the two schemes occurs when g 2 < 0. As shown in Fig. 
9 below, while the propagator in the former remains analytic (for g 2 > —1), singularity is 
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present in the latter. Thus, when employing the self-consistent RG prescription to analyze 
the fixed-point solution, it is desirable to expand about some $o where U^(^o) > 0. 




Figure 9. P£(z) as a function of z with a=ln2 and 6=2,3,4 and 10 for gi =— 0.461533<0. The inverse 
propagator becomes singular when z<^/— . 



APPENDIX C: ASYMPTOTIC EXPANSION WITH SMOOTH SMEARING 

In Sec. IV we see that the momentum integrations for the fixed-point solution in 
the self-consistent RG are divergent due to the negative fixed-point value g 2 . Since g 2 ~ 
—0.4615 is a "small" negative number, we inquire how good the approximation is to make 
an asymptotic expansion about g\ = 0. 

To explore the asymptotic limit, we first expand £/*($) about the origin, and the (3 
functions read 

(3 2 = -2g 2 - gJi 

04 = ~e<74 + 3^4/2 - g 6 Ii (C.l) 
(3 6 = (2d - 6)g 6 - 30^4/3 + 15<?4#6^2 - g%h, • • ■ 



where 



f 1 z d (t) _ f 1 dzz d - l {abz b e- az ) _ f 1 dzz d ~ l (,dp kb 

n[92) ~L [z*{t)+g 2 } n ~ Jo (z*+g 2 ) n ''Jo {^W 

(C.2) 

At M = 2 the 2x2 linearized RG matrix becomes 

f-2 + gtI 2 (gZ) -h(g* 2 ) \ 
M b =\ , (C.3) 

V -&gTh{gl) -e + SgthmJ 

where the non-trivial Wilson-Fisher fixed points are located by solving the transcendental 
equations 

n fi -*rf-*U Tf-^ (\+ z d (t)[(Q + e)g* 2 +ez 2 (t)} 

= 6g 2 I 2 {g 2 ) + eh{g 2 ) = / dt L — — -, (C.4) 

Jo (z 2 (t)+g* 2 ) 
22 



and g% = e/3J2(^|)- Note that in the limit b — > oo, I n — > (l+^2) _n , and we readily recover 
the sharp cutoff results. 

Expanding I n {g 2 ) about g 2 = 0: 



i"n(^2) = Oft / T- 

^0 1 



d-l+b-2n „-az b 



z + n — 1\ d-2(i+n) /b + d — 2(i + n) 



i=0 



n — 1 



)■ 



(C.5) 



we have for e = 1, 

(92> at) 



5ci - y/25cf + 44c c 2 6ci - yggcj + 44c c 2 \ 
22^ ' 3(cf - 4c c 2 ) / 



(C.6) 



and 



3(5ci-44c c 2 -cly / 25cf+44c^') 



T2T^f 



2c 2 (-6ci+ v /25c^+44coc 2 ) S 
3( C 2-4coc 2 ) 2 



(C.7) 



where c = a" 1 / b r(l + 1/6), d = a 1 / f T(l — 1/6), and c 2 = r(l-3/6). The b dependence of 
the exponent v is depicted in Fig. 10 below. From the Figure, we see that the expansion 
is only reliable for sufficiently large b. Indeed as b — > oo, v approaches 0.525 which is the 
sharp cutoff result for M = 2. 
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Figure 10. Asymptotic value of v as a function of b. 



The asymptotic limit with p~k,m = i+(pjh) m can ^ e ex pl° re d in a similar manner. In 
this case we have 

dzz d ~ x /,d,3 fc , m \ f°° dzz d + m - x 



j (— \ / azz ( u opk,m\ / azz"^" 

I n (92) = - Y^[ k ^J~) =m T~9 ^' C - 8 

Jo (z 2 +92) v dk J J (z 2 + g 2 ) (l + z m ) 2 
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and the fixed-point solution for M = 2 is found by solving 

. dz z m+2( 7 g* + z 2j 



Making an asymptotic expansion, we have 



m 



[l + z m ) 2 (z 2 + gl) 2 



(C.9) 



rt-\ V"V . l fi + n-l\n(d-2(i + n)) /n(d - 2(i + n)) 
In(g2) = 2^(-92) l [ _ )- = -csc( 



i=0 



77,-1 



m 



777 



(CIO) 



provided that d — 2(i + n) + m > 0. The resulting linearized RG matrix takes on the same 
form as Eq. (C.7), where Co = c\ = ^csc(^) and c 2 = 3^-csc(— ). The m dependence of 
v is shown in Fig. 11. Once more, the results are only reliable in the large m limit, where 
v w 0.525, in agreement with that obtained using the sharp cutoff. 
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Figure 11. Critical exponent v as a function of m. 
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